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We present a first measurement of CP asymmetries in neutral B decays to D + D~ , and updated 
CP asymmetry measurements in decays to D* + D~ and D*~ D + . We use fully-reconstructed decays 
collected in a data sample of (232±3) x 10 6 T(AS) -» BB events in the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. We determine the time-dependent asymmetry parameters 
to be S D , +D - = -0.54 ±0.35 ±0.07, C D ,+ D - = 0.09 ± 0.25 ± 0.06, S D ,- D + = -0.29 ± 0.33 ± 0.07, 
C D — D + = 0.17 ± 0.24 ± 0.04, S D+D - = -0.29 ± 0.63 ± 0.06, and C D+D - = 0.11 ± 0.35 ± 0.06, 
where in each case the first error is statistical and the second error is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charge-parity (CP) violation is described in the Stan- 
dard Model (SM) by a single irreducible complex phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) quark mix- 
ing matrix V Q]. The B-meson system provides an ex- 
cellent probe for testing the completeness of the CKM 
mechanism in a variety of CP asymmetries Mea- 
surements of CP violation in B° — > (cc)K ^ decays Q 
by the BABAR and Belle collaborations have 
precisely determined the parameter sin.2/3, where (3 is 
arg[-V cd y c b/V td V t ^]. The current world average of 
sin2/3 = 0.726 ± 0.037 is in good agreement with the 
range implied by other measurements in the context of 
the SM [6j , providing evidence that the CKM mechanism 
is the main source of CP violation in the quark sector. 

Decays of B° mesons to pairs of charged Z^*' mesons 
can also be used to determine sin2/3. These decays pro- 
ceed to leading order via a tree-level color-allowed b — > 
ccd transition. The presence of a gluonic penguin contri- 
bution with a different weak phase is expected to change 
the magnitude of the CP asymmetry by not more than a 
few percent 0). However, additional contributions from 
non-SM processes may lead to shifts as large as A/3 0.6 
in some models Q . Interference between SM penguin and 
tree amplitudes can additionally provide some sensitivity 
to the angle 7 = arg [ -V ud V:JV cd V; h } Q . 

In this Letter we present a first measurement of CP 
asymmetries in the decay B° — > D + D~ and improved 
measurements of CP asymmetries in B° — ► D* + D~ and 
B° -> D*-D+ decays |H Ej. The results are based 
on an analysis of (232 ± 3) x 10 6 T(AS) -> BB decays 
recorded by the BABAR detector 12j at the PEP-II e+e" 
collider. 

The selection of B° — > D* + D~ candidates is sim- 
ilar to that of our previous analysis ^3- We recon- 
struct D* + in its decay to D°ir + , where the D° is re- 
constructed in one of four final states: K~tt + , K~ir + ir°, 
K~ir + ir~ir + , or K®ir + Tr~~ . The D~ is reconstructed in 
the final states K + ir~ir~ or K®tt~. The K® candidates 
are reconstructed from 7r + 7r~ pairs within 15 MeV/c 2 of 
the nominal K® mass |_13j . The tt° candidates are re- 
constructed as photon pairs with an invariant mass be- 
tween 115 and 150 MeV/c 2 ; each photon must have en- 



ergy above 30 MeV in the laboratory frame and the sum 
of the photon energies must exceed 200 MeV. We re- 
quire the D° and candidates to have reconstructed 
invariant masses within 20 MeV/c 2 of their respective 
nominal masses, except for D° decays with a 7r° daugh- 
ter, which must be within 35 MeV/c 2 of the nominal D° 
mass. The B° — > D + D~ candidates are reconstructed 
solely through the decay of D T — > K ± tt t tt t . Charged 
kaons are required to be incompatible with a pion hy- 
pothesis on the basis of measurements of particle identi- 
fication quantities. 

To reduce background from continuum events 
(e + e~ — > qq, q = u,d, s, c) is reduced, we exploit the con- 
trast between the spherical topology of BB events and 
the more jet-like nature of continuum events. We require 
the ratio of the second-to-zeroth order Fox- Wolfram mo- 
ments 0| to be less than 0.6. We also use a Fisher 
discriminant, constructed as an optimized linear combi- 
nation of 11 event shape variables |l5|: the momentum 
flow in nine concentric cones around the thrust axis of 
the reconstructed B° candidate, the angle between that 
thrust axis and the beam axis, and the angle between 
the line-of-flight of the B° candidate and the beam axis. 
The Fisher discriminant selection requirement increases 
the signal significance by 2% in the case of B° — > £)* ± D =F 
and 9% in the case of B° -> D+D~. 

For each candidate, we construct a likelihood vari- 
able £ ma ss from the differences between the reconstructed 
masses and the nominal masses of the D* + , D + , and 
D° candidates 10]. The £ maS s variable is the prod- 
uct of the likelihood functions for the three candidate 
types. The likelihood for D + and D° is parametrized 
with a single Gaussian function, while the mass differ- 
ence rri£)*+ — m,£,o is parameterized as the sum of two 
Gaussian functions. The computed value of £ maS s and 
the difference AE between measured energy of the B° 
candidate in the center-of-mass frame and half the center- 
of-mass energy, AE = Eg — (y/s/2), are used to reduce 
the combinatoric background. Maximum allowed values 
for both — ln£ mass and \AE\ are set for each individ- 
ual final state separately, optimized using a Monte Carlo 
simulation 
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FIG. 1: Distribution of m ES for (a) W -> D* + D~, (b) 'fi' -> D*~D + and (c) "b' -> D + D~ candidates. The shaded areas 
represent the contributions from background events. The dashed and solid curves describing the background and signal plus 
background distributions respectively are explained in the text. 



nificance. 

The technique for measuring the CP asymmetries is 
analogous to previous BABAR measurements described in 
detail elsewhere 0] . After the reconstruction of a B° — > 
£)(*)±_D=F candidate B C p, we assign the remaining tracks 
in the event to the other B meson -Btag- We compute a 
proper time difference At and its estimated uncertainty 
a At from the reconstructed decay vertices of Bcp and 
Stag- The tracks assigned to -Btag are used to determine 
the Btag flavor and thus the flavor of the Bcp meson 
at At = [lflj . Events are classified in one of six tag 
categories and must have an estimated probability w of 
assigning the wrong flavor to B tag less than 45%. 

Taking into account the uncertainty in the vertex po- 
sition and tag flavor, the observed At distribution for 
50 d (*)± d t s ig na i events Fg p (Ai) is described by: 



F± ( At) 



.-MA 

4t b o 



1 ± (1 - 2w)[S f sin(Am d At') 



-C f cos(Am d At')] | ® i2(At - At'; a At ), (1) 

where the difference between the observed and true de- 
cay time differences At — At' is described by the em- 
pirical resolution function R(At — At';aAt)- This func- 
tion is parametrized as the sum of three Gaussians, a 
'core' and a 'tail' Gaussian, each with a width and mean 
proportional to a At, and an outlier Gaussian centered 
at zero with a width of 8 ps. The values of the B° 
lifetime t b o and the B°-B° oscillation frequency Am^ 
are fixed to (1.536 ± 0.014) ps and (0.502 ± 0.007) ps" 1 
respectively 0. We determine Sf and Cf sepa- 
rately for D+D-, D*+D-, and D*~ D+ . If only tree- 
graph contributions are present, we expect Sd+d- = 
— sin2/3; C D + D - = 0, and C D *+ D - = —C D — D +. Ad- 
ditionally, under these conditions we have S^+d- = 



-Xsin(2/3 
X = 



1 



r 2 



S) and S D *- D + = —Xsin(2(3 — S), with 
and where 5 is the difference of the 



-d+ 



strong phases for B° -> D*+ D~ and B° -> D*~ D+ . If 
the magnitudes of the amplitudes for B° — * D* + D~ and 
B° -> D*-D+ are equal Q, then C D , +D - = C D ,- D + = 
0. To determine the values of id for each of the tag 
categories and to increase the precision on the reso- 
lution function parameters, we simultaneously fit to a 
large sample -Bfl av of reconstructed B° decays to the 
flavor eigenstates h + (h + = n + ,p + , and af) and 

J/t/jK*°(K*° -> K+tt-) 

The beam-energy substituted mass toes = [( s /2 + 
Pi ■ Pb) 2 I 'Ef — Pq] 1 / 2 , where the initial total e + e~ four- 
momentum (Ei,pi) and the B momentum p*B are de- 
fined in the laboratory frame, is used to determine the 
composition of the reconstructed D^^D^ samples. We 
use only the region ibes > 5.2 GeV/c 2 , which includes 
a large sideband of pure background events. These 
events are included in order to determine the proper- 
ties of combinatoric background present in the signal re- 
gion. Backgrounds are incorporated with empirical de- 
scriptions of their At spectra. The backgrounds include 
prompt decays (associated with background from con- 
tinuum events), and non-prompt decays with a At de- 
scription similar to Eq. Both components are con- 
volved with a resolution function distinct from that of 
the signal, parametrized as the sum of two Gaussians. 
Based on Monte Carlo studies we expect a significant 
flavor asymmetry in the non-prompt background of the 
B° — > -D* ± D =F samples, because the _D* ± candidate is 
usually a true D*^ while the D is more often incorrectly 
reconstructed. This flavor asymmetry is parametrized 
via values of Cf and Sf of the non-prompt background 
that are allowed to vary in the fit. 

The At and Toes distributions are fit simultaneously. 
The toes distribution, shown in Fig. ^ allows a deter- 
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mination of a signal probability for each event. In signal 
events, the values of 771es accumulate near the nominal 
B° mass with a resolution of approximately 2.6 MeV/c 2 . 
The fitted ttjes shapes consist of a Gaussian distribution 
for the signal and an ARGUS function [2(| for the com- 
binatoric background. The total number of selected can- 
didates iV can d and the signal yield N s i g are shown in Ta- 
ble 0] From detailed Monte Carlo simulations of generic 
B decays, we expect some background events to peak in 
the m-ES signal region due to cross-feed from other decay 
modes. The fraction of events in the signal Gaussian due 
to this peaking background is estimated to be (7.0±6.2)% 
for B° -> D^D* and (13.6 ± 6.2)% for B° -> D+D~. 



CP asymmetry 



Sample 



AT, 



cand 



iVsig purity 



y° -f d*+d- 



y° ^ D+D - 



993 126 ± 16 0.49 ± 0.03 
1038 145 ± 16 0.49 ± 0.03 
538 54 ±11 0.37 ±0.06 



TABLE I: Candidates, signal yield and purity for each of the 
samples. The purity is defined as the fraction of signal events 
A/aig/ATcand in the region ?ties > 5.27 GeV/c 2 . 

The increase in statistics since our last measure- 
ment 10 for B a —>■ D* ± D T has allowed some refine- 
ments in the analysis. These include an improved treat- 
ment of signal probabilities as determined from the ities 
spectrum, and additional floating parameters for the de- 
scription of the background of the CP sample. We have 
also improved the event reconstruction, candidate selec- 
tion, and tag-flavor determination. The present effective 
tagging efficiency Q — 30.5% 0, a relative increase of 
5% over the algorithm previously used. 

We perform separate fits for each of the three CP sam- 
ples. There are in total 54 floating parameters describing 
the At distributions. These are Cf and Sf for signal (2) 
and background (2) , the average mistag fractions Wi and 
the differences Awi between B° and B° mistag fractions 
for each tag category i (12), parameters for the signal At 
resolution (7), parameters for background At distribu- 
tion (4) and resolution (3) of the Bflav and CP samples, 
and values for Wi and Awt for the prompt (12) and non- 
prompt (12) background of the -BfW sample. 

The likelihood fits yield the following results: 

S d *+d- = -0.54±0.35(stat.)±0.07(syst.), 

C D , +D - = 0.09±0.25(stat.)±0.06(syst.), 

S d —d+ = -0.29±0.33(stat.)±0.07(syst.), 

C D *- D + = 0.17±0.24(stat.)±0.04(syst.), 

S D+D - = -0.29±0.63(stat.)±0.06(syst.), 

C D+D - = 0.11±0.35(stat.)±0.06(syst.). 

Projections of the fit onto At for the three different 
CP samples are shown in Fig. |3 together with the raw 



a raw 
A CP 



(At) 



N + (At) ~ N^(At) 
N + (At) + N_{At)' 



(2) 



where N+(At) (N-(At)) is the number of B a -> 
£)(*)±£>t events with a B° (if) tag. 

The systematic uncertainties on 5/ and C / are sepa- 
rately evaluated for each of the decay modes. The domi- 
nant systematic uncertainty is the precision to which we 
are able to ascertain, using a Monte Carlo simulation, 
that the measurement method is unbiased (giving sys- 
tematic uncertainties in the range 0.03-0.06). Other im- 
portant uncertainties are due to the amount of peaking 
background and its potential CP asymmetry (0.01-0.02); 
assumptions on the At resolution function (0.01-0.03); 
and potential differences between the mistag fractions for 
the -Bfw and Bqp samples (0.01-0.02). Further sources 
of systematic uncertainty include the shape of the mES 
distribution, detector misalignment, uncertainty in the 
beam energies, and the possible interference between the 
suppressed b — > ucd amplitude with the favored b — > cud 
amplitude for some tag-side decays j^H. The total sys- 
tematic uncertainty is considerably smaller than in our 
previous measurement (0.10-0.14), primarily due to fewer 
assumptions about the background of the CP sample. 

In summary, we have performed a first measurement 
of CP asymmetries in the decay B° — > D + D~ . We 
have also updated our CP asymmetry measurements in 
B° — ► D* + D~ and B° — * D*~D + , superseding our pre- 
viously published results [lfj. The measured values are 
consistent with St — — sin2/3 and Cf = 0, expected in 
the SM for a tree-level-dominated transition with equal 
rates for B° -► D*+D~ and B° -> D*~ D+ . Since the 
dominant uncertainties are statistical, we anticipate im- 
proved precision with data collected in the future. 
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